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ABSTRACT 

We examine the constraints on high-redshift star formation, ultraviolet and X- 
ray pre-ionization, and the epoch of reionization at redshift z r , inferred from the 
recent WMAP-5 measurement, r e = 0.084 ±0.016, of the electron-scattering optical 
depth of the cosmic microwave background (CMB). Half of this scattering can be 
accounted for by the optical depth, r e = 0.04-0.05, of a fully ionized intergalactic 
medium (IGM) at z < zqp ~ 6-7, consistent with Gunn-Peterson absorption in 
neutral hydrogen. The required additional optical depth, Ar e = 0.03 ± 0.02 at z > 
^GP) constrains the ionizing contributions of "first light" sources. WMAP-5 also 
measured a significant increase in small-scale power, which lowers the required ef- 
ficiency of star formation and ionization from mini-halos. Early massive stars (UV 
radiation) and black holes (X-rays) can produce a partially ionized IGM, adding 
to the residual electrons left from incomplete recombination. Inaccuracies in com- 
puting the ionization history, x e (z), and degeneracies in cosmological parameters 
(O m , fife, o~8, n s) add systematic uncertainty to the measurement and modeling of 
T e . From the additional optical depth from sources at z > zqp, we limit the star- 
formation efficiency, the rate of ionizing photon production for Pop III and Pop II 
stars, and the photon escape fraction, using standard histories of baryon collapse, 
minihalo star formation, and black-hole X-ray preionization. 

Subject headings: cosmology: theory — cosmic microwave background — early 
universe — intergalactic medium 
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1. INTRODUCTION 

In recent years, an enormous amount of exciting cosmological data have appeared, accom- 
panied by theoretical inferences about early galaxy formation and the first massive stars. Many 
of these inferences were reactions to first-year (WMAP-1) results (Kogut et al. 2003; Spergel et 
al. 2003) from the Wilkinson Microwave Anisotropy Probe (WMAP). WMAP-1 inferred a high 
optical depth to the cosmic microwave background (CMB) and suggested early reionization of 
the intergalactic medium (IGM). Other conclusions came from simplified models for the stellar 
initial mass function (IMF), atomic/molecular physics, radiative processes, and prescriptions 
for star formation rates and escape of photoionizing radiation from protogalaxies. 

The CMB optical depth and other cosmological parameters have been refined significantly 
in the recent WMAP-5 data (Hinshaw et al. 2008). In this paper, we use these new measure- 
ments to constrain the efficiency of first-light ionizing sources. We focus on the reionization 
epoch, defined as the redshift z r when the IGM becomes nearly fully ionized over most of its 
volume (Gnedin 2000, 2004). Our knowledge about reionization comes primarily from three 
types of observations: hydrogen Lya absorption in the IGM, high- z Lya-emitting galaxies, and 
CMB optical depth. Optical spectroscopic studies of the "Gunn-Peterson" (Lya) absorption 
toward high-redshift quasars and galaxies imply that H I reionization occurred not far beyond 
z GP ~ 6 (Becker et al. 2001; Fan et al. 2002, 2006). Ultraviolet spectra suggest that He II 
reionization occurred at z ~ 3 (Kriss et al. 2001; Shull et al. 2004; Zheng et al. 2004). The 
detection of high-redshift (Lya-emitting) galaxies (Hu & Cowie 2006) suggests a somewhat 
higher redshift, z r > 6.5. Data from WMAP, after three years (Spergel et al. 2007) and five 
years (Hinshaw et al. 2008), suggest that reionization might occur at z ~ 10, with sizeable 
uncertainties in measuring and modeling the CMB optical depth. 

There appears to be a discrepancy between the two epochs, zqp pa 6-7 and z r pa 10. 
However, the WMAP and Lya absorption results are not necessarily inconsistent, since they 
probe small amounts of ionized and neutral gas, respectively. Both the H I absorbers and 
ionized filaments in the "cosmic web" (Cen & Ostriker 1999) are highly structured at redshifts 
z < 10 and affect the optical depths in Lya. In order to effectively absorb all the Lya radiation 
at z ^ 6 requires a volume-averaged neutral fraction of just xm ~ 4 x 10~ 4 (Fan et al. 2006). 
Simulations of the reionization process (Gnedin 2004; Gnedin & Fan 2006) show that the 
transition from neutral to ionized is extended in time between z = 5 — 10. The first stage 
(pre-overlap) involves the development and expansion of the first isolated ionizing sources. 
The second stage marks the overlap of the ionization fronts and the disappearance of the last 
vestiges of low-density neutral gas. Finally, in the post-overlap stage, the remaining high- 
density gas is photoionized. The final drop in neutral fraction and increase in photon mean 
free path occur quite rapidly, over Az pa 0.3. 
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Initial interpretations of the high WMAP-1 optical depth (r e = 0.17 ± 0.03) were based 
on models of "sudden reionization" , neglecting scattering from a partially ionized IGM at 
z > z r . The three- year (WMAP-3) and five-year (WMAP-5) data resulted in considerable 
downward revisions, with recent measurements giving r e = 0.084 ±0.016 (Hinshaw et al. 2008; 
Komatsu et al. 2008). These analyses assumed a single step to complete ionization (x e = 1 at 
z < z r ), while Dunkley et al. (2008) also explored a two-step ionization history with a broader 
distribution of r e . The WMAP-5 data alone gave a 5<r detection, r e = 0.087 ± 0.017 (Dunkley 
et al. 2008), while WMAP-5 combined with distance measurements from Type la supernovae 
(SNe) and baryon acoustic oscillations (BAO) gave r e = 0.084 ± 0.016 (Hinshaw et al. 2008; 
Komatsu et al. 2008). The single-step analyses gave a reionization epoch z r = 10.8 ± 1.4 at 
68% confidence level (C.L.). The x 2 curves for two-step ionization provide little constraint on 
the redshift of reionization or on an IGM with low partial ionization fractions, x e <C 0.1, as 
we discuss in § 3. Figure 3 of Spergel et al. (2007) and Figure 8 of Dunkley et al. (2008) show 
that the parameters x e and z r are somewhat degenerate, each with a long tail in the likelihood 
curves. This emphasizes the importance of contributions to r e from ionizing UV photons at 
z > z GP from early massive stars and X-rays from accreting black holes (Venkatesan, Giroux, 
& Shull 2001, hereafter VGS01; Ricotti & Ostriker 2004; Begelman, Volonteri, & Rees 2006). 

Models of the extended recombination epoch (Seager, Sasselov, & Scott 2000) predict a 
partially ionized medium at high redshifts, owing to residual electrons left from incomplete 
recombination. Residual electrons produce additional scattering, Ar e ps 0.06 between redshifts 
z ~ 10-700. These effects are computed in CMB radiation transfer codes such as CMBFAST 
and RECFAST, but only a portion of this scattering affects the large angular scales (£ < 10) 
where WMAP detects a polarization signal. X-ray preionization can also produce CMB optical 
depths r e > 0.01 (VGS01; Ricotti, Ostriker, & Gnedin 2005). 

As we discuss in § 2.1, a fully ionized IGM from z = back to z G p = 6.1 ± 0.15, the 
reionization epoch inferred from Gunn- Peterson absorption (Gnedin & Fan 2006), produces 
optical depth, r e m 0.041 ± 0.003. This represents nearly half of the WMAP-5 measurement, 
r e = 0.084 ± 0.016. If the epoch of Gunn-Peterson reionization is zqp = 7 (Wyithe et al. 
2008), the contribution is r e = 0.050 ± 0.003. Therefore, the high-redshift ionizing sources 
are limited to producing an additional optical depth, Ar e < 0.03 ± 0.02. In § 3, we discuss 
the resulting constraints on the amount of star formation and X-ray activity at z > 7. We 
also parameterize the efficiency of various star formation environments, as determined by halo 
mass and the metallicity of stellar populations, to address several questions. Can reionization 
occur at relatively low efficiency, corresponding to star formation in massive Lya-cooled halos 
(virial temperature T vir ps 10 4 K), or must it take place in H 2 -cooled minihalos (T vir Pi 10 3 
K)? Does one require the high ionizing efficiencies of Population III (metal- free) stars, or can 
reionization be achieved by lower-efficiency Population II stars? 
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The IGM is thought to have a complex reionization history, with periods of extended 
reionization for H I and He II (Venkatesan, Tumlinson, & Shull 2003; Cen 2003; Wyithe & 
Loeb 2003; Hui & Haiman 2003; Benson et al. 2006). Because r e measures the integrated 
column density of electrons, there are many possible scenarios consistent with the WMAP-5 
data. In addition, reionization models are quite sensitive to the amount of "small-scale power" 
for ionizing sources, which depend on cr 8 , the normalized amplitude of fluctuations. 



2. OPTICAL DEPTH TO ELECTRON SCATTERING 



2.1. Analytic Calculation of Optical Depth 



To elucidate the dependence of CMB optical depth on the epoch of reionization (redshift 
z r ), we integrate the electron scattering optical depth, T e (z r ), for a homogeneous, fully ionized 
medium out to z r . For instantaneous, complete ionization at redshift z r , we calculate r e as the 
integral of n e aTd£, the electron density times the Thomson cross section along proper length, 



T e (z r ) 



n e a T (l + z) 1 [c/H{z)\ dz . 



(i; 



We adopt a standard ACDM cosmology, in which (dl/dz) = c(dt/dz) = (1 + z)~ x \cj H(z)], 
where H(z) = H [Q m (l + z) 3 + ^a] 1 / 2 and Vt m + Q A — 1 (no curvature). The densities of 
hydrogen, helium, and electrons are written tin = [(1 — Y)Vt b p cr / run}^ + z) 3 , % e = ynn, and 
n e = nn{l + y), if helium is singly ionized. We assume a primordial helium mass fraction, 
Y = 0.2477 ± 0.0029 (Peimbert et al. 2007) and define y = (Y/4)/(i - Y) w 0.0823 as the He 
fraction by number. Helium contributes 8% to r e , assuming single ionization (He II) at z > 3. 
An additional r e 0.002 comes from helium reionization to He III at z < 3 (Shull et al. 2004). 
The critical density is p cr = (1.8785 x 10~ 29 g cnT 3 )/i 2 where h = (iZo/100 km s" 1 Mpc" 1 ). 
The above integral can be done analytically: 



T e (z r ) 



c 



2Q b 
3fi 

n 



p cr (l-Y)(l + y)a T 



[{Q m (l + ^) 3 + Q A }V2_ 1 ] 



= (0.00435)/i 70 [{0.2794/if 2 (l + z r f + 0.7206} 1/2 - l] , (2) 
where we write (1 — Y)(l + y) = (1 — 3Y/4) and use the updated WMAP-5 parameters, 



tt h h 2 



L b n = 0.02265 ± 0.00059 and Vt m h 2 = 0.1369 ± 0.0037. To this formula we add the extra 
scattering (r e 0.0020) from He III at z < 3. With this He III contribution, equation (2) 
yields r e = 0.040, 0.045, and 0.050 for z r = 6.0, 6.5, and 7.0, respectively. 



For large redshifts, Q m (l + z) 3 3> Cl A , and the integral simplifies to 



T e (z r ) ^ ( — 



2Q b 

3^1 m 



3Y/A)a T 



(1 + Zr) 



3/2 



(0.0522) 



1 + z r 



3/2 



(3) 
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Here, we have scaled to a reionization epoch z r ~ 7. From the approximate expression in 
equation (3), we see that r e oc (pcAn^Hu 1 ). Thus, T e is nearly independent of the Hubble 
constant, since p cr oc h 2 while the combined parameters, Q^h 2 and Q m h 2 , are inferred from 
D/H, CMB, and large-scale galaxy motions. The scaling with h therefore cancels to lowest 
order. A slight dependence remains from the small Q\ term in equation (2). 

If we invert the approximate equation (3), we can estimate the redshift of primary (Gunn- 
Peterson) reionization, (1 + z GP ) « (7.77)[r e (^ GP )/0.05] 2 / 3 , scaled to the value, r e = 0.05, 
expected for full ionization back to z GP 7. As discussed in § 2.2, this is approximately the 
WMAP-5 value of optical depth, r e = 0.084 ± 0.016, reduced by Ar e « 0.03. This additional 
scattering, Ar e , may arise from high- ,2 star formation, X-ray preionization, and residual elec- 
trons left after incomplete recombination. The latter electrons are computed to have fractional 
ionization x e m (0.5 — 3.0) x 10~ 3 between z = 10-700 (Seager et al. 2000). Inaccuracies in 
computing their contribution therefore add systematic uncertainty to the CMB-derived value 
of r e . Partial ionization may also arise from the first stars (Venkatesan, Tumlinson, & Shull 
2003, hereafter VTS03) and from penetrating X-rays produced by early black holes (VGS01; 
Ricotti & Ostriker 2004, 2005). For complete sudden reionization, Komatsu et al. (2008) esti- 
mated z r pa 10.8 ± 1.4 at 68% C.L., by combining WMAP-5 data with other distance measures 
(SNe, BAO). The WMAP-5 data alone (Dunkley et al. 2008) imply r e = 0.087 ± 0.017, with 
z r = 11.0 ± 1.4 (68% C.L.). Their likelihood curves allow a range 7.5 < z r < 13.0 at 95% C.L. 
They also claim that WMAP-5 data exclude z r = 6 at more than 99.9% C.L. 

The additional ionization sources at z > zqp will contribute electron scattering that 
may bring the WMAP and Gunn-Peterson results into agreement for the epoch of complete 
reionization. In our calculations, described in § 3, we make several key assumptions. First, 
we assume a fully ionized IGM out to z GP ps 6-7, accounting for both H + and ionized helium. 
Second, we investigate the effects of IGM partial ionization at z > z GP . Finally, in computing 
the contribution of residual electrons at high redshifts, we adopt the concordance parameters 
from the WMAP-5 data set. The CMB optical depth is formally a 5er result, which may 
improve as WMAP refines its estimates of the matter density, Q m , and the parameters, er 8 and 
n s , that govern small-scale power. Both <r 8 and n s have well-known degeneracies with r e in 
CMB parameter extraction (Spergel et al. 2007; Dunkley et al. 2008). The derived parameters 
may change in future CMB data analyses, as the constraints on r e continue to evolve. In 
addition, inaccuracies in the incomplete recombination epoch and residual ionization history, 
x e (z), add uncertainties to the CMB radiative transfer, the damping of £- modes, and the 
polarization signal used to derive an overall r e . 
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2.2. Residual Electrons in the IGM 



We now discuss the contribution of residual electrons in the IGM following the recombina- 
tion epoch at z ~ 1000. Scattering from these electrons is significant and is normally accounted 
for in CMB transport codes such as CMBFAST (Seljak & Zaldarriaga 1996) through the post- 
recombination IGM ionization history, x e (z). However, a number of past papers are vague 
on how the ionization history is treated and on precisely which electrons contribute to total 
optical depth r e . This has led to confusion in how much residual optical depth and damping 
of CMB power has been subtracted from the CMB signal. It is important to be clear on the 
definition of the integrated r e when using it to constrain the amount of high- z ionization from 
first-light sources. Modern calculations of how the IGM became neutral have been done by 
Seager et al. (2000), although their code (RECFAST) continues to be modified to deal with 
subtle effects of the recombination epoch and the atomic physics of hydrogen (ns — > Is and 
nd — > Is) two-photon transitions (Chluba & Sunyaev 2007). 

To illustrate the potential effects of high-z residual electrons, we have used numbers 
from Figure 2 of Seager et al. (2000), the top-panel model, which assumed a cosmology with 
Htot = 1, ^ = 0.05, h = 0.5, Y = 0.24, and T CM b = 2.728 K. At low redshifts, z « z r , just 
before reionization, they find a residual electron fraction x e $ ~ 10 -3 3 . We fitted their curve 
for log x e out to 500 to the formula: 



where a ps 2.303 x 10~ 3 . More recent recombination calculations (W. Y. Wong & D. Scott, 
private communication) using WMAP parameters (fi 6 = 0.04, h = 0.73, VL m = 0.24, f2 A = 0.76) 
and Y = 0.244 find somewhat lower values, x e $ ps 10~ 3 - 67 with a ps 2.12 x 10" 3 . The 
lower x e fi to the faster recombination rates arises from their higher assumed baryon density, 
Q b h 2 = 0.0213, compared to Q b h 2 = 0.0125 in Seager et al. (2000). 

To compute the electron-scattering of the CMB from these "frozen-out" electrons, we use 
the same integrated optical depth formula (equation 1), in the high-z limit, where (di/dz) = 
(1 + zy^c/Hiz)} ^ (c/# )fim /2 (l + z)~ 5/2 - We integrate over the residual-electron history, 
from z r ps 7 back to a final redshift Zf ^> z r , to find 



x e (z) = x efi lO - 001 ^*) ps (5 x 10- 4 ) exp[a(l + z)\ , 



(4) 



(Ar e ) 



res 



( _c_\ \ PcrQ- ~ Y)a T VL b 




(l + z) 1/2 x e (z) dz 



r (i+z f ) 



(3.18 x 10- 3 ) x e:0 \ 



J(l + Z r ) 



u 1 ' 2 exp(au) du . 



(5) 



A rough estimate to the residual scattering comes from setting a = and adopting a constant 



7 



ionized fraction x e (z), 



(Ar e ) 



res 



(2.12 X 1(T 3 ) ( Xe ) [(1 + Zff 12 -(1 + Z r 



) 3/2 ] • 



(6) 



This estimate gives r e = 0.039 for z r — 6, Zf — 700, and (x e ) w 10~ 3 . More precise values of 
r e can be derived from the exact integral (eq. 5) by expanding the exponential as a sum and 
adopting the limit Zf ^> z r , 



For z > 500, the approximate formula (eq. 4) underestimates x e , but one can integrate 
the appropriate curves (Seager et al. 2000) using piecewise-continuous linear fits. Table 1 
lists the fitting parameters, x 6i0 and a, for various redshift ranges, together with the extra 
contribution, Ar e . These calculations give a total optical depth in residual electrons Ar e xs 0.06 
back to z — 700. These electrons have maximum influence on angular scales with harmonic 
£ max w 2z 1/2 ks 20-50 (Zaldarriaga 1997). At higher redshifts, x e rises to 10 -2 - 1 at z = 800 
and to 10 -1 ' 1 at z — 1000, where the CMB source function will affect the "free-streaming" 
assumption used in CMBFAST (Seljak & Zaldarriaga 1996). 



The WMAP-5 measurements (Hinshaw et al. 2008; Komatsu et al. 2008) of fluctuations 
in temperature (T) and polarization (E) have been interpreted to estimate total electron- 
scattering optical depth of r e pa 0.084 ± 0.016. The central value comes from computing the 
likelihood function for the six-parameter fit to WMAP-5 data (TT, TE, EE) marginalized 
with BAO and SNe distance meaasures. Because of the challenges in translating a single 
parameter (r e ) into a reionization history, x e (z), it is important to recognize the sizable error 
bars on r e . At 68% C.L. (Figure 1 of Komatsu et al. 2008; Figure 6 of Dunkley et al. 2008), 
when marginalized against other parameters such as "tilt" (n s ), the optical depth ranges from 
r e = 0.06 - 0.11. The lower value is only slightly above the optical depth, r e = 0.05, for a fully 
ionized IGM back to redshift zqp ~ 7. The higher value, r e = 0.11, clearly requires additional 
ionizing sources at z > zqp- 

In § 2.1, we showed that ~ 50% of this r e can be accounted for by a fully ionized IGM 
at z < zqp- Observations of Lja (Gunn-Peterson) absorption toward 19 quasars between 
5.7 < z < 6.4 (Fan et al. 2006) are consistent with a reionization epoch z r ~ zqp = 6.1 ± 0.15 
(Gnedin & Fan 2006) . According to equation (2), this produces r e = 0.041 ±0.003, where our 
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error propagation includes relative uncertainties in z r (2.5%), Vl m h 2 (8.4%), and Vl b h 2 (4.0%). 
Detections of Lya emitting galaxies at z ~ 6.6 (Hu et al. 2002; Kodaira et al. 2003; Hu & 
Cowie 2006) suggest that the epoch of full reionization might be as high as z r ~ 7. This 
reionization epoch corresponds to r e m 0.05 and is consistent with the small neutral fraction to 
which Gunn-Peterson test is sensitive, particularly when one accounts for density bias in the 
observed high- ,2 ionization zones around quasars (Wyithe, Bolton, & Haehnelt 2008). Residual 
post-recombination electrons produce a substantial optical depth from z 10 back to z xs 700, 
which uniformly damps all angular scales. However, their effect on the TE and EE power is 
considerably less on large angular scales (£ < 10). Thus, we can identify a significant portion 
of the WMAP-5 observed optical depth through known sources of ionization. The "visible 
ionized universe" out to zqp = 6-7 accounts for r e = 0.04-0.05, while high-z partial ionization 
could contribute additional Ar e 0.01-0.03. For this study, we investigate the requirements 
to produce an additional optical depth, Ar e = 0.03 ±0.02, from star formation and early black 
hole accretion at z > zqp- 

Our calculations represent an important change in the derivation of z r from r e , and sug- 
gests that the amount and efficiency of high- ,2 star formation need to be suppressed. This 
suggestion is ironic, since WMAP-1 data initially found a high r e = 0.17 ± 0.04 (Spergel et al. 

2003) implying a surprisingly large redshift for early reionization, ranging from 11 < z r < 30 
at 95% C.L. (Kogut et al. 2003). These results precipitated many investigations of star for- 
mation at z = 10 — 30, some of which invoked anomalous mass functions, very massive stars 
(VMS, with M > 140 M ), and an increased ionizing efficiency from zero-metallicity stars 
(VTS03; Wyithe & Loeb 2003; Cen 2003; Ciardi, Ferrara, & White 2003; Sokasian et al. 2003, 

2004) . Tumlinson, Venkatesan, & Shull (TVS04) disputed the hypothesis that the first stars 
had to be VMS. They showed that an IMF dominated by 10 - 100 M© stars can produce the 
same ionizing photon budget as VMS, generate CMB optical depths of 0.09-0.14, and still 
be consistent with nucleosynthetic evidence from extremely metal-poor halo stars (Umeda & 
Nomoto 2003; Tumlinson 2006; Venkatesan 2006). 

Although the IGM recombination history, x e (z), is included in calculations of CMBFAST 
and in CMB parameter estimation, the best-fit values of r e from WMAP-5 and earlier CMB 
experiments have been attributed exclusively to the contribution from the first stars and/or 
black holes at z < 20. The contributions from post-recombination electrons (20 < z < 1100) 
have not always been subtracted from the data. This post-recombination contribution was 
relatively small in some earlier models of reionization (Zaldarriaga 1997; Tegmark & Silk 
1995) that explored optical depths of r e = 0.5 — 1.0 and suggested reionization epochs up 
to z r ~ 100. However, with current data indicating late reionization, it becomes particularly 
important to consider contributions to r e at z > Zqp, prior to the first identified sources of light. 
The latest WMAP-5 results (Hinshaw et al. 2008) find a lower r e , but they also suggest more 
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small-scale power available for reionizing sources, owing to higher normalization parameters, 
cr 8 ~ 0.817 ± 0.026 and Vt m h 2 ~ 0.1369 ± 0.0037. This change marks a significant increase 
(A<7 8 = +0.056) over WMAP-3 data. Between WMAP-1 and WMAP-3, the decrease in a$ 
was offset by a reduction in spectrum tilt. However, this index did not change significantly in 
the recent data, going from n s = 0.958 ± 0.016 (WMAP-3) to n s = 0.960^;^ (WMAP-5). 
Alvarez et al. (2006) argued, from the lower values of r e and a 8 , that both WMAP-3 and 
WMAP-1 data require similar (high) stellar ionizing efficiencies. Haiman & Bryan (2006) use 
the lower r e to suggest that massive star formation was suppressed in minihalos. Our results 
on a lower Ar e make these requirements even more stringent, as we now quantify. 

In the next two sub-sections, we consider two scenarios for producing partial IGM ioniza- 
tion at z > Zqp. The first considers ionizing UV photons from the first massive stars, which 
will create small "bubbles" of fully ionized gas inside ionization fronts. The second scenario 
examines extended partial ionization produced by X-rays from accreting black holes. 

3.1. Ionization by Hot Stars 

Semi-analytic and numerical models of reionization (Ricotti, Gnedin, & Shull 2002a,b; 
VTS03, Haiman & Holder 2003) show that the efficiency of ionizing photon injection into 
the IGM can be parameterized by the production rate of ionizing (UV) photons from mas- 
sive stars. In our approximate models, we compute the hydrogen ionization fraction as 
x e (z) = €\jy fb(z) / cl(z) , where euv = -Vy/*/csc and fb{z) is the fraction of baryons in col- 
lapsed halos. To account for the nonlinearity of recombinations, we define cl = (^hii)/ ( n Hii) 2 , 
the space-averaged clumping factor of ionized hydrogen. We assume that cl is the same for H II 
and He III and that N 1 f*f esc is constant with redshift. As discussed by Gnedin & Ostriker 
(1997), particularly in their § 3.1, clumping of the photoionized gas sharpens the reioniza- 
tion transition, with a recombination time much less than the local Hubble time at z ~ 6-8. 
However, the integrated optical depth through this epoch is not greatly affected by clumping. 

The main effects of clumping arise when we translate the production rate of ionizing pho- 
tons into the ionization history, x e (z). When the local density of photoionized hydrogen is 
close to ionization equilibrium, (^HiTm) ~ (®ip(T) ^e^mi), and the quadratic dependence of 
recombination rates reduces the local ionization fraction by a factor c^ 1 . When the medium is 
out of equilibrium, either ionizing or recombining, this approximation, x e (z) = euv fb( z ) I c l( z ) , 
is less accurate. To assess the validity of our approximation, we have run a set of reionization 
models, integrating differential equations for the propagation of ionization fronts (VTS03) to 
derive x e (z). When we apply our approximation with cl = 10, it closely approximates the 
integrated hydrogen ionization history with a clumping factor of Cl = 40. These two curves 
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are virtually indistinguishable for nearly all the redshifts approaching reionization. Thus, an 
artificial boost by a factor of 4 in the clumping factor would correct for any deviations in 
our approximation from the true ionization history. Our approximation is not dramatically 
different in shape from the true ionization history, at least from our semi-analytic calculation, 
and requires at most a constant multiplicative correction in the clumping factor. We are cur- 
rently pursuing numerical experiments to understand the relationship between clumping and 
reionization, using adaptive mesh refinement simulations (Hallman et al. 2007) and exploring 
different weighting schemes for clumping on the sub-grids. 

The free parameters in this model are encapsulated in the "triple product", f^N^f^, 
where /* represents the star-formation efficiency (the fraction of a halo's baryons that go into 
stars), N 7 is the number of ionizing photons produced per baryon of star formation, and / csc 
is the fraction of these ionizing photons that escape from the halo into the IGM. We can now 
use our calculations to constrain the amount of high-z star formation through the product 
of these three parameters. We solve for the ionization history, using our previous formalism 
(VTS03) in which the baryon collapse history, fb{z), is computed through the Press-Schechter 
formalism with the cosmological parameters from WMAP-5. The propagation of ionization 
fronts is followed through the production rate of ionizing photons minus recombinations. 

Each of these three parameters has some dependence on the halo mass and environment 
(Haiman & Bryan 2006; Ricotti & Shull 2000). Since we have already parameterized the 
intrahalo recombinations through / csc , we account for the loss of ionizing photons on IGM 
scales through Cl in two forms: (1) a power-law form with slope (3 = —2 from the semi- 
analytic work of Haiman & Bryan (2006); and (2) the numerical simulations of Kohler, Gnedin 
& Hamilton (2007), using their case C (overdensity 5 ~ 1 for the large-scale IGM) for Cr, the 
recombination clumping factor corresponding to our definition of c^. In the latter case, the 
clumping factor is almost constant (c^ ~ 6) until the very end of reionization. Together, these 
two different cases provide bounds on the range of possible values, although they are quite 
similar over the epoch z ~ 6-8 that dominates the scattering. In our calculations, we use the 
Haiman & Bryan (2006) formulation of Cl(z). 

With these assumptions, we can use equation (5) and the allowed additional optical depth, 
Ar e < 0.03, to constrain the ionizing efficiency of the first stars. In Figure 1, we plot curves of 
efficiency factor, euv> as a function of the resulting Ar e . The two panels illustrate the changes 
in required efficiency produced by the WMAP-5 increase in "small-scale power" compared to 
WMAP-3 parameters. These differences arise primarily from the amplitude of fluctuations, erg, 
which increased from a$ = 0.761+°;°^ (WMAP-3) to <j 8 = 0.817 ± 0.026 (WMAP-5). In each 
panel, we show two curves, corresponding to star formation in Lya-cooled halos, with virial 
temperature T vir 10 4 K, and in H 2 -cooled minihalos (T vir 10 3 K). For these calculations, we 
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adopt the clumping factors from Haiman & Bryan (2006), a star-formation efficiency /* = 0.1, 
and a range of escape fractions / csc = 0.1-0.4, The two curves (red and blue) illustrate cases 
designated as Population II (metal-enriched) and Population III (zero-metal) stars. Both 
models adopt Salpeter initial mass functions (IMF), for which we find: (1) iV 7 = 60, 000 for a 
metal-free IMF (10-140 M ) that agrees with both CMB and nucleosynthetic data (TVS04); 
and (2) iV 7 = 4000 for a present-day IMF (1-100 M ). Note that N y = 34,000 is consistent 
with r e « 0.084 ± 0.016 (WMAP-5) if photons arise from zero-metal stars (Tumlinson 2006). 
Values of iV 7 were derived (TVS04) from the lifetime-integrated ionizing photon production 
from various stellar populations and IMFs and used as inputs in reionization models. 

Figure 1 shows how the required additional optical depth, Ar e = 0.03 ± 0.02, translates 
into the efficiency required of massive stars forming in high- ,2 halos. These plots are meant to be 
used as indicators of the types of star-forming halos and the stellar populations that generate 
ionizing photons, per baryon that passes through stars. The two panels also illustrate the 
sensitivity of the efficiencies to cosmological parameters (WMAP-3 vs. WMAP-5), primarily 
the amplitude of fluctuations, cr 8 , which sets the small-scale power of halos at a given epoch. 
The interpretation of Figure 1 comes from the intersection of the two rising curves (red for 
massive halos, blue for mini-halos) with the two green bands, which mark the range of expected 
efficiencies of Pop II and Pop III stars. 

One can draw several general conclusions from Figure 1. First, the higher amplitude of 
small-scale power (larger erg) in the WMAP-5 results (right panel) shifts both red and blue 
curves to the right. Lower efficiencies of ionizing-photon production are required to produce a 
given additional optical depth, Ar e . Second, considering just the WMAP-5 parameters (right 
panel), we see that "massive halos", defined as those with T vir m 10 4 K, can produce optical 
depths ranging from Ar e = 0.005-0.018, for a range of Pop III efficiencies euv = 600-2000. If 
one invokes more plentiful mini-halos (T vir ss 10 3 K), reionization by Pop II stars can produce 
a similar Ar e = 0.004-0.016. These results suggest that partial reionization at z > z GP could 
be achieved by minihalos (T vir pa 10 3 K), with a mixture of Pop III and Pop II stars, or a 
transition from one to the other, for which the rising blue curve gives Ar e pa 0.02-0.06 for 
/esc = 0.1-0.4. For these minihalos, Pop II stars produce r < 0.02, but r pa 0.02-0.06 can be 
achieved either by metal-poor stellar populations in halos with unusually high /* or / esc , or by 
metal-free massive stars in halos with very low star formation rates and/or escape fraction of 
ionizing radiation. 

Another robust conclusion from Figure 1 is that mini-haloes (T vir = 10 3 K) with Pop III 
stars can easily account for the entire additional optical depth (Ar e ~ 0.03) needed to explain 
the WMAP-5 data. Even in cases where reionization begins with Pop III stars (efficiencies 
~ 10 3 ) and then transitions to Pop II (efficiencies ~ 10 2 ), the optical depths are significant. 
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This is an important astrophysical issue, since the duration of the epoch dominated by metal- 
free massive first stars is still uncertain (TVS04). These results mark a significant departure 
from similar curves for a WMAP-3 cosmology (left panel), owing to the increased small-scale 
power in a WMAP-5 cosmology relative to WMAP-3. Although the scalar spectral index n s 
did not change appreciably from year 3 to year 5 of WMAP, both the matter density, Q m , 
and the normalization, as, rose by several percent, an appreciable effect for power available 
for low-mass halos. Thus, reionization occurs slightly earlier in a WMAP-5 cosmology (Az r ks 
1-2) in our models, lowering the required efficiencies by a factor of about 4. 

Our predictions can be compared to those of Haiman & Bryan (2006), who used WMAP-3 
data to constrain the efficiency factor, euv, relative to a fiducial value e m j n i = 200 for minihalos. 
They argued that the efficiency for the production of ionizing photons must have been reduced 
by an order of magnitude, in order to avoid overproducing the optical depth. However, their 
constraints were based on the WMAP-3 value, r e xs 0.09, whereas half that scattering is 
accounted for from the fully ionized IGM at z < zqp- In our formalism (Figure lb), the 
efficiency constraints are less severe, and Pop III minihalos can easily produce Ar e 0.06 
with the expected efficiencies, euv ~ 10 3 . Much larger efficiencies, close to those of metal-free 
stellar populations, and large values of /* and / esc are required for larger halos to produce 
Ar e ~ 0.02. It makes some difference what form is assumed for the clumping factor, cl, in 
constraining the ionizing efficiency of the first stars. The level of shifts in the (red and blue) 
efficiency curves is comparable to the shifts between WMAP-3 and WMAP-5. 



3.2. Ionization by Accreting Black Holes 

The CMB optical depth also constrains the level of X-ray preionization from high-redshift 
black holes. Ricotti et al. (2005) were able to produce large optical depths, up to r e w 0.17, 
using accreting high-z black holes with substantial soft X-ray fluxes. Their three simulations 
(labeled M-PIS, M-SN1, M-SN2) produced hydrogen preionization fractions x e = 0.1 - 0.6 
between z = 15 and z = 10, with large co-moving rates of star formation, (0.001 — O.1)M 
Mpc -3 yr _1 , and baryon fractions, ujbu ~ 10~ 6 - 10 -5 accreted onto black holes. Any significant 
contribution to r e from X-ray preionization requires x e > 0.1. Therefore, the lower value 
of r e from WMAP-5 reduces the allowed X-ray preionization and black-hole accretion rates 
significantly compared to these models. Ricotti et al. (2005) find that the IGM at z > zqp was 
highly ionized (x e ^> 0.01). In this limit, there were few X-ray secondary electrons, and most 
of the X-ray energy went into heating the ionized medium. By contrast, the WMAP-5 optical 
depth suggests that the IGM was much less ionized at z > zqp- 

Relating the effects of X-ray ionization from early black holes to a Ar e and a related X-ray 
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production efficiency factor is less straightforward compared to the star-formation case, for the 
following reasons. First, unlike ionization by UV photons, X-ray ionization is non-equilibrium 
in nature, and the timescales for X-ray photoionization at any epoch prior to z = 6 typically 
exceeds the Hubble time at those epochs (VGS01). Therefore, it is more difficult to establish 
a one-to-one correspondence between X-ray production at an epoch and the average efficiency 
of halos. In addition, X-ray ionization (whether from stars or black holes), at least initially 
when x e < 0.1, will be dominated by secondary ionizations from X-ray- ionized helium electrons 
(VGS01) rather than from direct photoionization. This may constrain the physical conditions 
in the IGM (e.g., the level of He ionization) rather than those in the parent halo, when we 
attempt to translate a Ar e into an X-ray ionization efficiency. Thus, it may be difficult to 
make precise inferences about the black hole density and accretion history from r e . 

We can define an X-ray ionization rate and efficiency parameter, analogous to the previous 
case for massive star formation. We define x e (z) = exfb(z), with no clumping factor, where 
ex = NxfBnfesc and fb(z) is the fraction of baryons in collapsed halos. Here, /bh is the 
average fraction of baryons in black holes in halos at z > 7 and Nx is the number of X-ray 
ionizations of the IGM produced by X-rays and secondary ionizations, per baryon accreted 
onto such black holes. As we now describe, /bh may approach 10 -4 , and standard accretion 
estimates give Nx ~ 10 6 photons/baryon, so that ex ~ 100. 

Ricotti & Ostriker (2004) and Ricotti et al. (2005) define a parameter o;bh = 10~ 7 to 
10~ 6 , equal to the fraction of baryons that go into seed black holes. Current observations 
(e.g., Haring & Rix 2004) suggest that the BH-to-bulge mass fraction in modern galaxies is 
M B H/-^"buige ~ 1-4 x 10~ 3 . Thus, if the bulge-to-halo mass ratio is 0.065, /bh ~ 10~ 4 , after BHs 
have grown from the initial fraction through accretion, thereby producing X-rays ionization at 
high-2. To estimate N x , we assume Eddington-limited accretion at 10% efficiency and adopt 
a mean photon energy ~ 1 keV: 

(12) ~ (10 6 photons/baryon) E'^ . (8) 

Here, we assume that each X-ray photon produces ~ 12 hydrogen ionizations, primarily 
through secondary ionizations from X-ray photoelectrons (Shull & van Steenberg 1985; VGS01). 
In the partially ionized IGM, the free electrons come from H + , He + , and trace He +2 . We as- 
sume that the clumping factor (cl) and escape fraction (/ eS c) are roughly unity for X-rays, 
given their high penetrating power relative to UV photons. Figure 2 shows the allowed ad- 
ditional optical depth, analogous to the constraints of Figure 1, for X-ray efficiency in both 
Lya-cooled halos and minihalos. Evidently, X-rays from black holes located in high-redshift 
minihalos can produce Ar e ps 0.02, with a substantial ionizing efficiency, ex ~ 100, that rivals 
that of Pop II star formation. 



0.1m p c 2 
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In summary, we have shown that the revised ( WMAP-5) values of CMB optical depth, 
r e = 0.084 ± 0.016, lead to a more constrained picture of early reionization of the IGM. 
Approximately half of the observed r e comes from a fully ionized IGM back to zqp ~ 6- 
7. The additional Ar e at z > Zqp probably arises from the first massive stars and from 
accretion onto early black holes. Some of the observed r e may come from scattering from 
residual electrons left from recombination; inaccuracies in computing this ionization history 
add systematic uncertainty to the CMB inferred signal. We have assumed extra scattering, 
At 6 = 0.03 ± 0.02 at z > z GP and used this to constrain the efficiencies for production and 
escape of ionizing photons, either by ultraviolet photons from the first massive stars (Fig. 1) 
or by X-rays from accretion onto early black holes (Fig. 2). 

In both cases, the picture is of a partially ionized IGM at redshifts z = 6 — 20. For X-ray 
pre-ionization by early black holes, equation (6) can be used to provide an estimate of the 
effects of partial reionization. Between redshifts z 2 ~ 20 and z\ 6, an IGM with ionized 
fraction x e $ = 0.1 (Ricotti et al. 2005) would produce (Ar e ) = (0.018)(x ei o/0.1), which is a 
significant contribution to the observed r e = 0.084 ± 0.016. Ricotti et al. (2005) suggested a 
large x e ~ 0.1 — 0.6 and pushed their black-hole space densities and accretion rates to large 
values in order to reach the WMAP-1 estimates of r e = 0.17. Because such large values of r e 
are no longer required, the black hole densities and IGM ionization fractions are likely to be 
considerably less. 

All these constraints depend heavily on uncertain parameterizations of the efficiency of 
star formation and ionizing photon production. The most sensitive of these parameters is cr$, 
although the details of the clumping factor are comparably important. With more precise 
measurements of CMB optical depth from future missions, there is hope that more stringent 
constraints on high-z star formation and black-hole accretion will be possible. These data 
include additional years of WMAP observations, as well as the Planck mission. 
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Table 1. Ionization Fractions and Residual Optical Depth 



Redshift Range 


x e,0 a 


a a 


Ar e b 


7 < z < 500 


5.00 x 10~ 4 


2.303 x 10~ 3 


0.0256 


500 < z < 600 


2.94 x 10~ 4 


3.224 x 10~ 3 


0.0135 


600 < z < 700 


1.28 x 10~ 4 


4.606 x 10~ 3 


0.021 



a Parameters for x e (z) = x ej oexp[— a(l + 2;)] (see text) 

b Residual integrated optical depth computed over 
given redshift range (see eq. 5) 
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Fig. 1. — Locus of efficiency factors, euv = -^7/* /esc, for production and escape of photoion- 
izing radiation vs. additional optical depth, Ar e , from hot stars at z > zqp- Models assume 
cosmological parameters from WMAP-3 (left panel) and WMAP-5 (right panel). WMAP-5 
found more small-scale power (higher as) for mini-halo ionization, hence lower required effi- 
ciencies. Efficiency, euv = f*N 1 f esc , is defined (§ 3.1) for star formation in Lya-cooled massive 
halos (red curve, T vir > 10 4 K) and in H 2 -cooled minihalos (blue curve, T vir > 10 3 K). Green 
horizontal bands correspond to expected efficiencies for a Salpeter stellar IMF: euv ~ 600- 
2000 for metal-free Pop III (10-140 M ) and e uv « 40-150 for metal-enriched Pop II (1-100 
M & ). We adopt star-formation efficiency /* = 0.1 and photon escape fraction / esc ranging from 
0.1-0.4 (see § 3.1 for details). The intersection of the rising red curve (right panel) with the 
upper green band shows that Pop III halos (T vir « 10 4 K) can produce Ar e rs 0.01-0.02. The 
intersection of the rising blue curve with the two bands shows that mini-halos {T W1T « 10 3 K) 
can produce Ar e = 0.01-0.02 (Pop II) and Ar e up to 0.06 (Pop Ill/Pop II). 
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Fig. 2. — X-ray production efficiency, e x = NxfBnfesc, required to produce additional optical 
depth, Ar e , at z > zqp- We assume that / csc = 1 for X-rays. This efficiency is the product of 
the baryon fraction in black holes times the total ionizations produced per accreted baryon. 
Curves are for WMAP-3 parameters (left) and WMAP-5 (right). The required efficiencies 
are lower for WMAP-5 parameters, with their additional small-scale power. Two curves are 
shown, for halos with virial temperatures of 10 3 K and 10 4 K. Horizontal curves show a fiducial 
value, ex ~ 100, equivalent to /bh ~ 10~ 4 , N x ~ 10 6 (see § 3.2), and / esc = 1- 



